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Abstract 

In the case of polynomial potentials all solutions to 1-D Schrodinger equation 
are entire functions totally determined by loci of their roots and their behaviour 
at infinity. In this paper a description of the first of the two properties is given 
for fundamental solutions for the high complex energy limit when the energy 
is quantized or not. In particular due to the fact that the limit considered is 
semiclassical it is shown that loci of roots of fundamental solutions are collected 
of selected Stokes lines (called exceptional) specific for the solution considered 
and are distributed along these lines in a specific way. A stable asymptotic limit 
of loci of zeros of fundamental solutions on their exceptional Stokes lines has 
island forms and there are infintely many of such roots islands on exceptional 
Stokes lines escaping to infinity and a finite number of them on exceptional Stokes 
lines which connect pairs of turning points. The results obtained for asymptotic 
roots distributions of fundamental solutions in the semiclassical high (complex) 
energy limit are of a general nature for polynomial potentials. 
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1 Introduction 



As it is already well known [1, 2] the fundamental solutions (FS) [3] have appeared 
to play main role in one dimensional quantum mechanics (or in a multi dimensional one 
allowing a reduction to the one dimension) with analytic potentials, i.e. polynomial, 
meromorphic etc. In particular they allow us to solve all basic problems typical for 
the field - eigenvalue problems, scattering problems, problems of decaying and JWKB 
and adiabatic limits [4]. They are exceptional also among all the solutions to the 
corresponding Schrodinger equations because of their property of being Borel summable 
for the polynomial or meromorphic potentials [5]. The latter property allows us to 
recover these solutions from their JWKB series, both exactly and approximately if 
the semiclassical series is abbraviated in the latter case. Finally they allow us to 
use in effective calculations their approximations by JWKB formulae having these 
approximations under full control. 

Not surprising therefore that the fundamental solutions are also very important 
for mathematicians, used by them under the name of subdominant solutions, to study 
properties of solutions to Schrodinger equations (SE) for polynomial and meromorphic 
potentials [6]. Among many problems to be under considerations is the one of loci 
of zeros of solutions to SE as a function of coefficients of polynomial or meromorphic 
(rational) potentials (see for example [7]). This problem is also important from the 
physical point of view to mention the known relation between a number of real zeros 
of the quantized solution to SE and a number of an energy level corresponding to 
the solution. In a recent paper of Bender et al [8] a method of looking for zeros of 
eigenfunctions of eigenvalue problems with non-hermitian potentials was suggested as 
a tool for checking a possible completeness of the full set of such eigenfunctions. 

In another recent paper of Eremenko et al [9] the high energy limit has been con- 
sidered to study this problem for a quantized energy in polynomial potentials. The 
authors have shown that in this limit the problem simplifies greatly so that it can 
be standardized and the loci of zeros of the corresponding quantized FS's (i.e. these 
which are a solution to an eigenvalue problem) are exactly on Stokes lines, called ex- 
ceptional by the authors, of a Stokes graph (or global Stokes lines - a name used by 
mathematiciens) corresponding to a considered problem. 

Another semiclassical limit has been considered by Hezari [10] who investigated 
the problem of complex zeros of eigenfunctions of SE with real polynomial potentials 
of even degree in the limit H — > 0, where Ti is the Planck constant, while the energy 
parameter E was kept fixed. 

However, Hezari's eigenfunctions are as such for a problem of the quantized Planck 
constant (being by this a little bit unphysical). In fact the two cases, i.e. the energy 
quantized while the Planck constant is kept fixed and the energy kept fixed but the 
Planck constant is quantized have two different semiclassical limits for the quantized 
parameters, i.e. the high energy limit and the small fr-limit lead to different behaviour 
of the corresponding Stokes graphs and to different sets of eigenfunctions. Nevertheless 
since both these limits are of the same semiclassical nature at least mathematically it is 
therefore not surprising that Hezari's results of complex zeros eigenfunctions problem 
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are similar to those of Eremenko et al. We have discussed and generalized Hesari's 
results in a separate paper [11] 

In this paper we would like both to generalize the results of Eremenko et al [9] and 
to make the corresponding theorems more precise in the following aspects: 

1. to show that they are valid for not quantized FS's in the considered limit so that 
quantized cases can be seen as particular results of these general ones 

2. to show how the quantization procedure modifies unquantized zeros distributions 
of FS's 

3. to show that the exceptional Stokes lines of Eremenko et al [9] can be identified 
with the boundaries of the canonical domains corresponding to FS's 

4. to find the limit distributions of zeros of FS's along exceptional SL's using the 
explicite form of FS's and their high energy semiclassical limit. 

In fact our analysis provides a full description of the root distribution problem of 
FS's to SE with polynomial potentials in the high (complex) energy limit. 
The paper is organized as follows. 

In the next section the high energy limit of SE for polynomial potentials is consid- 
ered and it standard form is established. 

In section 3 Stokes graphs (SG) for the standardized high energy polynomial po- 
tential are considered and their possible standard forms are established in this limit. 

In the next section we define the full set of fundamental solutions for the standard 
high energy limit potential (— iaz) n and formulate the basic lemma on possible posisions 
of roots of the FS's. 

In section 5 two theorems are formulated establishing the precise positions of roots 
of FS's in the high energy limit. 

In section 6 the results of the previous section are extrapolated to an arbitrary 
polynomial potential rescaled correspondingly to the limit considered. 

In section 7 the quantized cases of FS's and an influence of the quantization on 
distributions of zeros of FS's are investigated. 

We conclude and summarize the results in section 8. 



2 High energy limit of Stokes graphs for polynomial 
potentials 

Consider the stationary 1-D Schrodinger equation with a polynomial potential 
P n (z) = a n z n + ... + aiz, a n ^ 0, n > 1, (a free term is assumed to be absorbed by the 
energy parameter E). We have: 

4>"{z) - |?(P B (*) - E)4>{z) = (1) 
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To standardize our problem we make the substitution z — > — ia (J^j n z, a — 1 for 
n odd and a — 1, e~ l « for n even so that we get for ip a (z) = ^— " -2^ 

^(z)-\ 2 W n (z,\)Mz)=0 (2) 

where 

W n (z, A) = (-iaz) n " 1 + E 6n-fc(-ia)^A-^(-iaz) n -* 

fc=i 

k 

2 2 2m n+2 a n -k f2m\^+2 

X=-a— T a n n E™, Qn-fc = ~ ^ti+2~ 1 — ) , fc = l,...,ra-l (3) 

ft 2(n+2) \ Tl J 

where the factor — ia has been introduced for a convenience what will be self-explaining 
later. 

Let us remind that for the case considered the Stokes graph (SG) is created as a set 
of all lines (Stokes lines (SL)) emerging from each root Zi(X), i — 1, ...,n, of W n (z, A) 
and satisfying one of the following equations 

m(\[ Z yJW n (£,\)dA =0, i = l,...,n (4) 

The roots ^j(A), % — 1, ...,n will be called also turning points. 

From (3) it is obvious that for \E\ —> oo, i.e. |A| — > oo, the limit Stokes graph 
corresponding to the rescaled problem is determind by the roots (turning points) of 
W nta (z) ee (— iaz) n — 1 and the phase of A. It is also clear that for the considered limit 
the roots 2* (A) are all simple and their loci are close to the ones of W nta (z) so that 
the SG corresponding to W n (z, A) differs only slightly from the one corresponding to 
W n>a {z) and can be obtained from the latter SG by its small deformation. 

By the above definition the Stokes lines and the corresponding Stokes graph are 
defined in the considered limit on the two sheeted Riemann surface i?2 with the turning 
points of W n!a (z) as the branch points of this surface. However since on these two sheets 

the values of \JW n)Cl (z) differ by a sign only the projections on the z-plane of the Stokes 
lines defined on each sheet coincide. 

Therefore considering a pattern of SL's on the cut z-plane C cut with cuts emerging 
from the turning points of W n ^ a (z) we see that the SL's on C cut are quasi continuous 
on the cuts despite the fact that they are pieces of different SL's collected from the two 
sheets of R 2 . 

In general, an SL emerging from a given turning point Zi(X) can run to infinity of 
C cut or end at another turning point Zj(X). An SL with the last property is called the 
inner one. 

An SG is called critical if at least one of its SL's is the inner one. It is called not 
critical in the opposite case. 
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3 Properties of Stokes graphs corresponding to the 
potentials (— iaz) n 

Consider now the SG's corresponding to W n ^ a (z). For this goal assume A to be 
real and positive for a while. 

For a = 1 there is a root of W n>a (z) at z = zq = i while the remaining ones are 

- 2&7T 

regularly distributed on the circle \z\ = 1 being located at z = Zk = ie l ~ , k = 



±1, ±2, ±[^^-] and at z = z* = — i for even n, so that the corresponding pairs of 
them satisfy the relation z k = —z~- k (where bar over z means its complex conjugation), 
i.e. these pairs are located symetrically with respect to the imaginary axis. 

For a = e~ l ™ all zeros are given by pairs Zk = e l ~ , k = ±1, ±2, ±|. 

It can be shown by a direct calculation (see Appendix 1 and [6]) that 

2i ain m B (l^) for a=l, fc=l,.., fc<" 



rz_ k , i -sin — d Hj, - for a — 1, ft = l, ft < tt 



where B(x,y) is the betha function. 
Therefore 
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I*" JwZXz)dz = ^ ft = 1,2..., (6) 



for any such a pair, i.e. each pair lies on the same Stokes line. 

Let us remind further that from each root of W n , a (z) (all the roots are simple) 
emerge three SL's. If these roots are members of a pair Zk,z-k, ft = 1,2..., then one 
of this lines runs from one root to another while the remaining pairs of SL's run to 
infinity of the z-plane. 

If a root of W nt0l (z) is not a member of any pair of them then all three SL's which 
emerge of it run to the infinity of the z-plane. 

Each pair of neighbour SL's emerging from the same root and running to the infinity 
forms a sector while the SL's alone lie on its boundary. 

However for an odd n there is still a sector which boundary is formed by the 
neighbouring SL's runnig to infinity and emerging from the last pair (with the highest 
value of ft = ! ^) of roots and by the SL linking this pair. 

It is easy to note that 2n + 4 is the total number of sectors lying on R 2 . However 
because of coincidence of SL's projected on C cut we can consider on this cut z-plane 
quasi sectors formed by these projected SL's. There are now n + 2 of such quasi 
sectors which will be called again sectors for simplicity. We can enumerate them 
correspondingly to roots attached to their boundaries. 

So for each pair of turning points Zk, z_k, k — l, 2..., the corresponding sectors are 
denoted by S k and S_k respectively while the single sector formed by the last pair of 
roots in the odd-n case is denoted by SWi and by Sn+2 in the even-n case and in the 
same case the single sector formed by the first pair of roots is denoted by So- 

If there are single roots, at z = z = i or z = z-± = —i, then the two sectors 
connected with the first root we denote by Sn+s (the left one) and S_n+3 for n odd and 



4 



by Sn+2 (the left one) and S_n+2 for n even. The second root at zn can exists only for 
n even and the two sectors connected with it are denoted by Sa (the left) and S-a. 




Fig. la SG for an odd n. Broken lines are cuts 

Now consider SG's corresponding to not standard cases of zeros distributions, i.e. 
when the standard distributions are rotated by some angle 0, < <p < -. Then 
the corresponding Stokes line linking a pair z^Zj, < argz| < argz' < 2ir, of not 
neighbouring roots of W ntOL {ze 1 ^) is determined by the equation 

m. f A J'J sJwZJ^)dz\ = (7) 

Making a rotation z — > ze tuJ where uj = |(arg z[ + arg z'j ± 7r) (we have to choose one 
of the possibilities) we come back to the standard configuration of roots of W n>a r(z) 
(perhaps with the different a') and the condition (7) reads then 

3? (Xe™ f'~ k JwZMdz) = (8) 

where Zk, Z-k is a pair of roots of W nj0l i{z) into which the pair z[, Zj is transformed after 
the cu-rotation. 

Taking now arg A = —uj we can satisfy (8) according to (6). 

We can conclude therefore that for an arbitrary not standard distribution of roots 
a not neighbouring pair of them can be connected by a Stokes line by a proper choice 
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of argument of A and if this choice is made we can cosider the arising Stokes graph in 
its standard position by making the corresponding ^-rotation. 

k Imz 




Fig. lb SG for an even n and a — 1. Broken lines are cuts 



We have got even more. Namely, if the condition (7) is satisfied for a pair z' i: z'j 
then it is satisfied also for all pairs of turning points which are symmetric with respect 
to the line argz = ^(arg^ + arg^) since the corresponding cu-rotations are the same 
for such pairs and coincide with the one for the pair z[, z'j. 

Consider now again the SG's corresponding to the standard configurations of roots 
and for real A. They are shown on Fig.la-e. What happens when A acquire a nonva- 
nishing phase (31 It is well known that in such a case each three SL's emerging from 
the same root rotate by an angle — ^ and of course for (3 = ±7r the whole SG comes 
back to its form for real A. 

However patterns of SG's arising by this /3-rotations are not all topologically differ- 
ent because of symmetric distributions of turning points. Namely, if we consider the SL 
linking a pair z*., (call it L^) and start to rotate the SG by sufficiently small f3 > 
then, since all SL's of the graph will rotate anticlockwise, Lj, will split into two SL's 
now runnig to the infinity each. Still enlarging (3 we can cause one of these splitting 
SL's emerging from z\. to meet the root Z-k-i, i>e. the two last roots find themselves on 
the same SL. But this is just the situation we have discussed previously so that by the 
corresponding tu-rotation in the z-plane we can return to the standard configuration 
of SG. It is easy to see that for this case the cj-rotation has to be done by one of the 
angles ±-, n ± -. A choice of the two latter angles corresponds to an odd n only and 
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has to be done only when really the SL arises between the roots z& and 

A Imz 




Fig.lc SG for an even n and a = 1 



Therefore as it follows from the above discussion it is enough to rotate the stan- 
dard SG by angles f3 chosen from the interwals (— -, -) to handle all topologically 
nonequivalent configurations of SG for the standard distributions of roots. 

Let us put therefore A = |A|e i/3 , < \(3\ < ^, in the condition 

&(ajT y/W n>a (S)df) =0 (9) 

where k takes all possible values following from the values of n, n = 2, 3, 4, ... . 

It then follows from the above discussion that all SG's corresponding to the A's 
chosen are non-critical, i.e. each of their SL's defined by (9) runs to the infinity of the 
z-plane, see Fig. 2. 

We shall show in the next sections that for the potentials {—iaz) n and W n (z) + 1 
in the limit |A| — > oo all zeros of each of their fundamental solutions are distributed 
along a boundary of its canonical domain, while this boundary is a collection of Stokes 
lines. 
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Fig. Id SG for an even n and a = e « 



4 Fundamental solutions for the potential (— iaz) n 

The fundamental solutions (FS) of the equation 

^j"(z)-\ 2 ((-iaz) n -l)^(z) = (10) 

are solutions defined on R 2 separately in each of 2n + 4 sectors. They are subdominant 
in the sectors which they are defined in, i.e. they vanish for \z\ — > oo inside the sectors. 
They can be given explicite forms of the following functional series [3]: 

where z € Sk and z& is a turning point lying on the boundary of Sk while 

W nj0l (z,z k )=[ Jw n , a (y)dy 

Jzk 

Xa ,k{z, A) = 1 + E (-^)V ai *;„(z, A) (12) 
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Fig.le SG for an even n and a = e" 



and 



Y 



a,k:n 



>,A) 



dyi / c?2/ 2 



hk(.Vn-l) 

_ e -2o-fcAVF n , a (z,j/i)^ ^ _ g-2cr fe AW'n lQ1 (l/l,10) 



dy n u> a (yi)u} a (y 2 ) . • • w«(2/ n ) x 



,-2(T| ! AB'n,a(!)»-l,Jn) 



n > 1 



with 



5 (*) 



16 
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(13) 



(14) 



The signatures er^ = ±1 present in the formulae (11)-(13) are defined in each 
particular sector Sk in such a way to ensure the inequality ^i(akXW nja (z, Zy.)) < to 
be satisfied in this sector. 

The integration paths 7fc(z) in (13) which start from the infinities of the correspond- 
ing sectors are chosen in such a way to satisfy the inequality $l(crkXW nta (yk, Vk+i)) > 
0, yk, yk+i £ 7, for each factor of the integrations in (13). Such paths are called canon- 
ical and a point z for which 7/t(z) can be chosen to be canonical is also called canonical 
with respect to the sector Sk- The latter property is completely of a topological nature. 

According to our earlier discussion we can assume the argument of A to vary in the 
interval (--.+-). 
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Let us note however that there is no necessity to consider the fundamental solutions 
on the whole R 2 - It is enough to define them and consider on C cut only because of the 
following reasons. 

As we have mentioned earlier there is one to one correspondence between the sectors 
lying on the C cut which can be considered as a one sheet of R 2 and the sectors lying 
on the second sheet C' cut of R 2 the latter sheet being then a complement of C cut to R 2 
connected with C cut by cuts. The correspondence between two such sectors is built by 
their coincidence when C' cut is projected on C cut . 

Let S C C cu t and S' C C' cut be a pair of such sectors. Then the two fundamental 
solutions defined in each of these two sectors coincide up to a constant. This coincidence 
is visible in the form given by (11)- (14) by the invariance of (JkW n ^ a {z, z k ) and a k u a {z) 

when passing from the sector S to S' since then W nj0t (z, z k ) and cu a (z) change their 

_ 1. 

signs simultanuously. Only the common factor W n ,a (z) of the solutions changes slightly 
under this operation acquiring one of the phase factors ±i. 

Therefore we can consider the fundamental solutions defined only on C cut and the 
solution ip a ,k{z, A) given by (11)-(14) is defined in the sector Sk the latter being enu- 
merated in the way described earlier. However if ip a} k(z,X) is defined in the sector 
crossed by a cut of C cut we have to remember about necessary changes in the forms 
(11)- (14) of the solution described above to keep its identity when the cut is crossed. 
This latter note is valid also when such a cut is crossed by any of the fundamental 
solutions if the latter is continued analytically along a canonical path. 

A collection Dk C C cut of all points canonical with respect to Sk is called a canonical 
domain corresponding to Sk- Of course Sk C Dk for each k [3]. 

A general rule for a given SL to belong to dD k is that a canonical path 7fc(z) when 
z approaches the line has to cross another SL emerging from the same turning point. 

For non-critical SG's (i.e. all SL's of which run to the infinities) dDk is composed 
of all turning points and of single SL's emerging from these turning points, see Fig.2. 

In other cases (i.e. for critical positions of SG) all three SL's emerging from a 
turning point can belong to dDk- The latter point is then just the root z_k joined by 
one of its Stokes lines with the turning point Zk, i.e. by the line which start from the 
former point and end at the latter. Such a line is then called an inner one. All points 
of the sector S-k cannot be then joined with the infinity of Sk by cannonical paths and 
therefore do not belong to D k . 

Each dD k is specific for its canonical domain D k . 

In what follows we will consider also a C cut (e)-plane arising from C cut by depriving 
the latter e- vicinities of the turning points Zk, k — 1, .., n. Namely, for any e > denote 
by Afc(e) the following circle vicinities of turning points z k : A fe (e) = {z : \z — z k \ < 
e}, fc=l,...,n. ThenC cut (e)=C cut \ULiA fe (e). 

Let L r k denotes a connected set of SL's contained in dD k and emerging from the 
turning points z r . Therefore for the non-critical SG's all L r k contain only single SL's and 
are disjoint pairways while for the critical SG's a unique difference with the previous 
case is connected with L\ and L k k which are not disjoint, the first one containing the 
inner SL between z k and z_ k while the second one containing all the three SL's emerging 
from Z-k, therefore also the inner SL between z k and Z-k- The lines L r k , r — 1, n, will 
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be called (after Eremenko et al [9]) exceptional with respect to the solution ij) a ,k{z)- Of 
course they are exceptional in that none of their points can be reached by the solution 
ipa,k{z) if the latter is to be continued to them along canonical paths. 

For a given U k let us denote by V k r (e) an e- vicinity of this set of SL's defined by the 
following conditions: 

1- Li C V k r (e) 

2. a boundary of V k (e) consists of (at most two) continues lines an Euclidean dis- 
tance of which to L r k is smaller than e (see Fig. 2 and 3.) 

Let us further denote by D ke a subset of D k given by D ke = D k \ V&(e) where 
It is clear that 14(e) is an e- vicinity of dDk (see Fig. 2 and 3.) 




Fig. 2 Vfc j£ - vicinity of a domain where ipk{z) cannot be continued canonically. 

The non-critical case 

Needless to say for the non-critical SG's every turning point zi belongs to V&(e) 
together with exactly one of its three SL's which emerge from it. In the case of the 
critical SG's the same is true for all z\ excluding z~k which belong to Vk(e) together 
with its three SL's if Zk is connected with Z-k by the inner SL. 

The following property of FS's is the key one for our further considerations. 
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Lemma 

In the domain Dk, e the factor Xa,k{z, A) of the solution (11) satisfies the following 
bound 

\Xa,k( z > A) - 1| < e~ - 1, |A| > A 
C a . e = liminf / |a; a (£)d£| < oo (15) 

7fe(je), zeDfc, e , k=l,...,nJ lk (z) 

where 7^(2) are canonical. 

We have left the proof of Lemma to Appendix 2. 




I 

Fig. 3 Vk, e - vicinity of a domain where cannot be continued canonically. 

The critical case 

There is a direct relation between above Lemma and so called semiclassical ex- 
pansions of x-factors of FS's (11) for A — > 00 but fixed z which will be needed in our 
further considerations. Such expansions have been considered in our earlier papers (see 
[4], ref.l and [5], ref.5) and have been given the following exponential forms: 
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where 

Y a , k ;o(z) = 1 
r z _i / —I fVp _i 

Ya,kA z )= dy n W n £{y p ) [W n £{y n ) dy^WnA (y p -i) x 

L^-|(y 2 ) £ dj/iw^i^o^i^i))"...) ] , p>i (17) 



and 



Note that X ajP (y), p > 1, are sector independent being point functions on the 
C cu t-plane which are given by the following recurrent formula (see [4], ref.2): 



X a;1 ( Z ) = U( Z ) = Wnj(z) (Wn}(z)J = W„j(z) 



n,a 



Z 



X«M = -lWn,l(z)W^ a (z)X a , p _ 1 (z) + 

Wn,l(z) *«,**«,m-k-l + X,p-l(*)) , P = 2, 3, ... (19) 

where C/ 2 n-2(-2 ; ) is a polynomial of the 2n — 2-degree. 

It follows from (19) that X a;2m , m > 1, have only poles at the turning points while 
2m+ii w > 0, have there the square root branch points. Therefore the same are 

the properties of Z+(z, A) and Z~(z, A) at these points respectively where Z£(z, A) + 

cr k Z-(z,X) = Z a , k (z, A) and 



/ 1 \ Zm 



2m 

2mUJ 



/1\ 2m+l 

Z-(*,A)=£ (— J (20) 



m>0 



If we now take into account that Xa,i-*jW = li m z->oo., Xa,^- 2 ) A) — Xa,j-n(ty where 
ooj is the infinite point of the sector Sj communicated canonically with the sector Si 
(see ref.l of [4] and ref.5 of [5]) then we get 

J^ { Zi + a t Z-) dz = J^{Z+ + * jZ - a)dz (2i) 

Since however o; L = —Uj then we get from (21) 

/ Z+(z,\)dz = (22) 
for any pair of canonically communicated sectors. 
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However the integration in (22) is now not limited by canonical paths since under 
the integral there are now no exponentials limiting this integration to canonical paths, 
i.e. these paths can be freely deformed with the integral beeing still convergent. It is 
easy to see that because of that a given integral (22) can be deformed to any other 
integral between any pair of sectors (i.e. not necessarily communicated canonically) as 
well as to any integral along an arbitrary loop. It means that residua of Z+(z,\) at 
the poles which it has at the turning points vanish. 

Therefore we conclude that the Riemann surface of Z+(z, A) is just the C-plane on 
which it is meromorphic with vanishing residua at its poles. It means of course that 
this is the property of each X Q)2m , m > 1 as well. Thus when integrating Z a ^(z,X) 
along contours starting and ending at the same points we get only contribution from 
the odd part of Z a;k (z, A), i.e. from o~ k Z~(z, A). 



5 High energy limit of loci of zeros of fundamental 
solutions for the potential (— iaz) n 

It follows from the Lemma that making A sufficiently large, i.e. A 3> C e , we 
can make Xa,k(z, A) arbitrarily close to unity for z G D k ^ and |A| > Ao- But this means 
that for such conditions ip a ,k(z,X) vanishes nowhere in D k)6 . Therefore, we have the 
following theorem about loci of zeros of i/) a ,k(z, A). 
Theorem 1 

For sufficiently large A all zeros of ip a)k (z, A) can lie exclusively in the completion 
C cut \ -Dfc, e of the domain D kje . 

In the case of a non-critical SG this completion coincides of course with V^ e while 
in the critical case it can contain also a whole sector. This is because in the case of 
the potential considered a connected set of SL's which contains more than three lines 
with three of them emerging from the same turning point z k has to contain also one 
more turning point, i.e. Z- k . If we consider a solution ip a ,k{z, A) defined in the sector 
Sk then the sector S- k and its boundary dS- k cannot be then connected with S k by 
canonical paths, i.e. they lie in the completion C cu t \ D k . We see that S- k is just the 
sector which has been mentioned earlier as contained in C cut \ D ke (see Fig.3). 

To formulate a theorem giving the precise positions of the roots of ip a ,k( z , A) men- 
tioned in Theorem 1 for the cases of non-critical SG's let us make the following several 
notes. 

The first one is that in choosing C cu t we can always choose the corresponding cuts in 
such a way to ensure all the SL's lying in C cut to be strictly (i.e. not quasi) continuous 
in this cut plane, i.e. none of these lines can cross any cut. With such a choice for 
every root z k of W njCe (z) the values of 3 (A f z Zk sjw n , a (y)dy} have definite sign on every 
SL emerging from z k depending on the line. 

The second one is, as we have already mentioned, that in the case of non-critical 
SG's every turning point z\ belongs to dD k together with exactly one of its three SL's 
which emerge from it and dD k is collected exactly of only such SL's. 
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The third one is the most important since it establishes the way of taking the limit 
A — > oo. Namely, as it follows from App.2 zeros of i> a ^{z.,X) scattered around such 
an exeptional line sliding down in the limit considered in directions to the respective 
turning point from which the line emerges. If these zeros are not sufficiently far from 
the turning point they unavoidably approach their limit just in this point. However 
in such a case the conditions of the validity of the representations (11)-(14) and their 
semiclassical expansions (16)-(19) are broken, i.e. the limit positions of zeros are outside 
of the domains D k>€ engaged in taking the limit. A necessary condition to get the 
limit loci of zeros different from the turning point is to take their initial loci to be 
A-dependent. This dependence in the formulae given below is not only necessary but 
also sufficient giving finite loci of zeros different from turning points. 

This dependence however defines also the way of getting the limit A — > oo inforcing 
to put |A| = [|A|] + A, < A < 1, where [|A|] is a step function of |A| (i.e. an integer 
not greater than |A| itself) and to consider the limit A — > oo conditioned by the fixed 
value of A. This kind of limit will be called regular to distinguish it from the free limit, 
i.e. with no conditions. Obviously, each A defines a different way of getting the limit 
A — > oo by ip a ,k( z , A) and its zeros as well. 

Nevertheless we can also consider zeros the limit loci of which are just turning 
points when A — > oo. However we have to remeber in such cases that the asymptotic 
formulae used to get these loci have to be bounded to the domains D k) e- 

The following theorems give the precise limit positions of roots of ip a ,k(z, |A|e J/3 ) for 
|A| — > oo up to all terms of A -1 . 

Theorem 2a 

Zeros CS(A), |A| = [|A|] + A, I = 0, ±1, ±2, q = 0, 1, 2, r = 0, ±1, ±2, .., 
°f ipa,k{ z , |A|e l/3 ), < \f3\ < -, i.e. in the non-critical cases, in the regular limit 
[|A|] — > oo ; are distributed on C cut uniquely along the corresponding exceptional SL's 
according to the formulae: 

/«!>,> 5^»<* A >) * = ± («i |A| l + r - i) X < 23 > 

where Ki{C^ r {\)) is a contour which starts and ends at C^qr(A) rounding the turning 
point z\ anticlockwise. 

Zeros (i k q r (X) for q > have the following semiclassical expansion: 

C(A) = E iC,(A) (24) 

p>0 A 

with two first terms given by: 

, -^JW n>a {y)dy = / l ' qr ° y/W n , a (y)dy = ±qme~ ip 



1 i-rc-W 



Ci(A) = Mr ~ ?A - t) . /Ak) (25) 

4 /m5 
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For q = we have instead C/,oro(^) = z i an< ^ 

r 2i +C i ( o ) T . 1 (A)/A ! i i n 

W n>a (y)dy= (r- -) — 



/ 



l A lr 
r > — lim sup 

71 \4>\<k 



zi+ee* 

W n , a (y)dy 



(26) 



as well as 



TTie signs ± a&ove are to fre chosen the same as of S (e l/3 \Jw^Jjj)dyj on the 
corresponding exceptional lines. 

First of all what follows from the formulae (24) and (25) is that, according to (4), 
C/orOJ 9 = ^'2) •••) ne on SL emerging from while the same can not be said about 



the first term of the expansion (24), i.e. the integral \ z l { qv yW n ^ a (y)dy is not purely 
imaginary if the first and next orders of A -1 are included in the expansion (24). It 
means that zeros Cz'gr(A) deviate from SL considered when higher order terms in (24) 
are included beginning from the first one. 

We get therefore the following picture of loci of zeros (^(A) in the considered regu- 
lar limit. For each A, < A < 1, these zeros are distributed along the SL emerging from 
Zi. For a given q, q = 1,2, every two neighbouring zeros are separated from each 

other by ine~ if3 / (XyJw n)0l (CkiSj ■ For a S iven 9, 1 = 1, 2, and r, r = 0, ±1, ±2, 
when A changes from zero to unity these zeros ocupy a line: 

which crosses the exceptional SL considered at z — C/qi-o ^ or ^ — ( r ~ \)/q- 

On the other hand for higher q — 1,2, a distance d q between the neighbouring 
such lines of zeros mesured along the exceptional SL (it is L\ this time) is equal to 

d q = me~ l13 1 y^ n ,a(C/,goo) an d vanishes with q — > oo. 

Fig. 4a illustrates the exceptional lines for ip a n+2 (z) for even n and a = e~ l ™ oc- 

-— i 

cupied by zeros of this solution with the distinguished exceptional line on which 

2 

the discribed above details of the zeros distribution are shown. 

To formulate the corresponding theorem for the critical cases of SG's, i.e. for real 
A, we have to notice that C cut \ = dD^ U S_k so that dDj, contains the point 
z^k together with their three SL's which emerge from it. This is just a set L^ k of 
exceptional SL's emerging from z^. 

Let us note further a role played in the critical cases of SG's by the integral 
A fz k ~ k \/Wn,a(y)dy = Aifc. For A real and positive it can always be represented as 
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IkX = — (r + R)m where r > is an integer and \R\ < |. Both r and R depend, 
of course, on A and due this dependence R can take on an arbitrary value from its 
domain of variation. Conversely, fixing R on some Ro we choose in this way an infinite 
sequence \k r (Ro) = — + Ro), r = 0,1,2,..., of A's growing to infinity and such 
that R(Xk, r (Ro)) = -Ro- We shall maintain further for such sequences a description 
regular as well as for the limit A — > oo itself using these sequences if \R\ < | while we 
call them singular if \R\ = |. As previously all other sequences of A's not restricted by 
any condition will be again called free. 

It should be stressed however that as it follows from the form of the x - f ac tors 
of FS's and what will be seen later when semiclassical limit A — > oo of particular 
formulae will be taken that to apply regular limits one needs to statisfy the inequality 
cos(i?7r) ;§> ^ > 0. However this condition is always satisfied for any \R\ < | when A 
is sufficiently large. 

Almz 




Fig. 4a SG for an even n and a = 1. The bold lines are exceptional SL's along which 
zeros of ^™+2 (z) are distributed in the regular limit A — >• oo. 
On one of them details of such ditribution are shown 

In Theorem 2a we have used regular sequences of A conditioned by fixed value of 
A. In Theorem 2b below the regular limit A — > oo will be also considered when R 
has a fixed value. This is essentially the main difference between the critical and non 
critical cases in taking the limit A — > oo. 

Theorem 2b 

Zeros CS(A), |A| = [|A|] + A, / = 0, ±1, ±2, I ^ -k, q = 1,2,.., r = 
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0, ±1, ±2, of ip a! k(z, \X\), i.e. in the critical cases, in the regular limit [|A|] — > oo ; 

1. e. with fixed A, are distributed on C cut uniquely along the corresponding exceptional 
SL 's according to the formulae: 

/«!>,> (fV^W - Tx z "^ A >) * = ± (« + r " i) X < 29 > 

where Ki{(j k q r {X)) is a contour which starts and ends at Ci# r W rounding the turning 
point Z[ anticlockwise. 

Zeros (i k g r (X) have the following semiclassical expansion: 

CS(A) = E iC(A) (30) 

p>0 A 

wif/j two first terms given by: 

J- rrr. , r, iH.qro 



J (k) ^W n , a {y)dy = J l ' qr ° \jW n , a {y)dy = ±qm 

= + (r -a\ - -) 

4 iSSS 



d(A) = ±(r - ?A - . S (31) 



For q = we have instead Cairo (A) = 



£ +<; " <A,/i ^) rf! /=( r -i)f 

+ee ^ 

\IW n , a {y)dy 



'Zl 

l A lr 
r > — hm sup 

7T UKtt 



(32) 



as well as 



8 y/WnAzi + tiSnW/X) 



The signs ± afrowe are to be chosen the same as of S ( ^W ntOC (y)dyj on the 
corresponding exceptional lines. 

In the case I = k the number q is bounded, i.e. q < \h\/^ and in (31) the minus 
sign is to be chosen according to our earlier conventions. 

Additionally in a regular limit Xk )S — > oo, i.e. with R fixed, where Afc )S = —^jj^m = 



[Xk,s\ + Afe )S (i2), s = 0,1,2,..., there are two infinite sequences of zeros C-l%^ 1 = 
1,2,..., r = 0, ±1, ±2, of i/) a ,k(z, \X\) distributed along the two infinite SL's of the 
sector S-k according to the following rules: 

/ k ' qr \/W nta {y)dy = - (q[X k ,s] + r -A + ^IY~ + aY~ f Za > kdy ± 

J- .„ 2 cos (ftr + ia Z^y) ± J- /.^ Z^y (34) 
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with the following first coefficients of the corresponding semiclassical expansion ofC-l% : 



W nia (y)dy = -qiTi 

1 R 1, /„ A 17T 



C ( -l%i(R) = - (r - ?A*,.(i2) - i + ¥ =F ^ ln2cos(ifor)) -=J 



(35) 



where the plus sign corresponds to a vicinity of the SL being the upper boundary of S_ k 
while the minus one to a vicinity of its lower boundary. 

Fig. 4b shows all the exceptional SL's of the odd n case occupied by zeros oitj} a ^{z, A) 
in the limits considered in the above theorem. 

All the above three theorems have been proved in Appendix 2. 

It follows therefore from the notes made above that the cases R = ±| have to be 
considered separately. It should not be surprising that they are just the quantized 
cases of the solution ip a , k (z, A) considered. 

6 High energy limit of loci of zeros of fundamental 
solutions for the potential W n {z, A) + 1 

2 

Let us denote — iaz by ( and A ™+ 2 by r] and let |argA| < ^. Then for the 
polynomial W n (z, A) we have: 



W n (z, A) -> W n (C, V ) = Wniux-^r,- 3 ?) = C ~ 1 + E K-kV k C~ k 

k=l 

2k 

b' n _ k = b n - k (-ia)—* (36) 

Let us note that bounding |argA| by ^ we always get the standard SG's corre- 
sponding to W n:a (z) as a limit of the corresponding SG's related to W n (z,\) when 
A — > oo. 

It follows from the above form of W n ((,i]) that if Co is its root for rj = then, by 
the implicite function theorem, there is also such a root ((rj) of W n ((, rj) which is close 
to Co f° r V sufficiently close to zero, i.e. £(77) = Co — ^n-i^ + ^K 7 ? 2 ) f° r V ~^ 0- ^ means 
of course that for sufficiently large A all the roots of W n (z, A) are simple and are close 
to the corresponding roots of W n ^ a (z), i.e to each root z k of the latter potential there 
is the root z k (X) of the former one such that z k (X) = z k — i-^b'^^^ + 0(\~~) 
for |A| — > oo, k — 0, ±1, .... 

Let us re-examin for the potential considered the respective notions of SG's and FS's 
and remaining ones related to them. A necessity to do it is caused by the A-dependence 
of the potential W n (z, A) + 1 and as it is seen from (36) this dependence is by powers of 

2 

A ™+ 2 . Such a dependence complicates an analysis of A-dependence of SG's, FS's and 
other notions. It is easier to make the respective analyses by a substitution A A^ in 



n-1 



19 



the corresponding Schrodinger equation (2) with the potential W n (z, A) + l by which the 
potential itself convert to depend on natural powers of A -1 and making simultanuously 
the same substitution for the SE with the potential W n , a (z) + 1. 



A Imz 




Fig. 4b SG for an odd n and a = 1. Bold lines are exceptional SL's along which 
zeros of ipn±2 (z) are distributed in the regular limit A — > oo 

The corresponding SG's can be now defined by using the new A-parameter for 
both the potentials. It does not disturb much properties of SG's for the potential 
W n>a {z) + 1 except that it makes these graphs more sensitive for A-changes (because 
A is now powered by (n + 2)/2), i.e. the graphs come back to their initial forms when 
A is rotated by the angle <fi = 2-n j(n + 2) rather than by it as previously, for example, 
and this is just this new angle which is now a period for such graphs' changes. 

In the case of the potential W n (z, A) + 1 the corresponding graphs defined by (4) 
where A is substituted by depend now on A and if A changes it needs a total 
change equal to 2tt of its argument to make the graphs coming back to their initial 
positions. Changing absolute value of A alone also forces the graphs to change while 
the graphs corresponding to W n>a {z) + 1 are insensitive on such a change of A. 

All these do not prevent us however to define the corresponding FS's by the formulae 
(11)-(14) together with the corresponding definitions of canonical paths, canonical do- 
mains, etc, including also constructions of domains Dk )E contained in canonical domains 
Dk together with e- vicinities Vk, e of their respective boundaries dDk, k = 0, ±1, and 
with exceptional SL's related to them. We have only to remember that all the above 
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notions depend now on A and can change with it continuously. 

Nevertheless despite the mentioned A-dependence Lemma of sec. 4 remains valid 
also in the case of the potential W n (z, A) + 1. This can be argued as follows. 

Consider a fundamental solution ip k (z, A) to the SE (2) corresponding to the poten- 
tial W n (z, A) + 1 and a respective domain D k ^ t (X) where this solution can be continued 
canonically to any of its point and let |A| > A with A sufficiently large to find each 
z k (X) inside corresponding circle vicinities A fc (e), k — 1, ...n. 

By the same arguments and with obvious changes the inequality (15) in the domain 
-Dfc,e(A) can be written for the solution ip k (z, A) as follows: 



\Xk(z,X) - 1| < exp 



C e (X) 

n+2 
\ 2 



C e (A)= liminf f \u(£, A)df| < oo 

7fc(«), zeCcut(e), fe=l,...,n+2 J*y k (z) 



(37) 



Therefore Lemma of sec. 4 sounds now (see App.2 for a proof): 
Lemma' 

In the domain D ke {X) the factor Xk( z ,X) of the FS solution tp k (z,X) constructed 
for the potential W n (z, A) + 1 satisfies the following bound 



\Xk(z, A) - 1| < exp 



a 



n + 2 

An 2 ' 



1, |A| > A (e) 
C e = C a . € + e 



(38) 



with A (e) sufficiently large and with fixed but arbitrary small e 

Therefore Theorem 1 remains valid also for the potential W n (z,X) + 1 with no 
changes. 

To formulate analogues of Theorem 2a and 2b we have to establish the corre- 
sponding forms of the semiclassical expansions (16)-(18) for the case of the potential 
W n (z, A) + 1. These expansions have to be a little bit different because of a dependence 
of the potential on A. 

First we note that the x-factors of the FS's ipk{z, A) with the potential W n (z, A) + 1 
and A changed as mentioned earlier satisfy the following equation: 



W n ~ 4 (z, A) ( W n Hz, X) Xk (z, A) ) + 2a k X^ X ' k (z, A) = 



(39) 



which can be transformed to the (pseudo-) integral equation of the form 

z,X) = l--^ [ Z WnHy,X)(WnHy^)Xk(z,X))"dy (40) 
2A 2 Joo k \ j 



Xk{z, 



corresponding to the Xfc-factor of the FS defined in the sector S k (X). 

An attempt to get a solution for Xk from (40) by iterations leads to a divergent 
series which however appears to be just a semiclassical expansion for Xk( z ,X) of the 
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form: 



X*(*,A)~xr(*>A) = E( 

p>0 v 



2A^ 



p _ 



(41) 



where Y kyP (z, A) are given by (17) with W n>a (z) substituted by W n (z, A). 

If further using (17) (with the substitutions mentioned) we make the following 
expansion: 



<?>o 



(42) 



then (41) can be given the following forms depending on a parity of n: 

n = 2m 



P>0 2^g=0 \{m + l)p+q ) 



J_ ^2m+2 Yk;2 P ,g(z) 



(43) 



1 ^2m+2 Yk;2p+l,q{z) 
~>2p l^q=0 ^(2m+3)p+ 9 : 

where the coefficients of the last expansions are given by the ones of the expansion (42) 



V -L- X^^m+l ik;2p+l,qK z ) „ _ <) | 1 

2m+3 Z^p>0 2 2 P <^<?=0 A( 2m+3 )P+9 ' ~~ 

^ 2A 2 



as: 



Yk;pq{z) — J2r=0 ( if) ^ / fe;r,(rrt+l)(p-r)+ ( ?(^) ) 
^ / fe;2p, ( 3(^) =Sr=o2 2( ' P r ^A;;2r,(2m+3)(p-r)+<?(^) 
n;2 P+ l, g (^)=E^o2 2(p - r) n ;2r+l,(2ra+3)(p-r)+<j(^) 



n = 2m 
n = 2m + 1 



(44) 



The exponential representations of x% s { z -> A) have similar forms to (16), i.e. 

X a k s (z,X)=exp( f Z Z k (y,X)dy 



Z k (z,X) = < 



n = 2m 



y^m ^fc;pq( z ) 

Z^p>l 2^<j=0 \(m+l)p+q 1 

E l -sp 2m + 2 Jy fc;2p,'?( z ) 
p>l 2 2p 2~iq=Q _\(2m+3)p+ 9 

O-fc J_ V^2m+2 ^fc;2p+l,q(z) _ 9 | I 

2m+3 Z^p>0 2 2 P ^q=0 \(2m+3)p+q 1 11 ~ Z '" t ' 1 

2A 2* 



> = 



X r(^,A) dz 



(45) 



However the corresponding recurrent relations for X k - pq (z) are now much more 
complicated. Nevertheless since Z k (z,X) satisfy the equations: 

Wn* (Wn^" + 2 (Wn* (Wn^j + d k X^ Z k + Wn* (Z\ + ^) = (46) 

then the partition Z k (z,X) = Z + (z,X) + a k Z~(z, A) still can be done with the same 
properties (22) for Z + (z, A) and the corresponding conlusions about its analytical prop- 
erties on the z-plane. 
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Unlike the case of the potential W na (z) + 1 we have to use also the following 



W n Jz) 



asymptotic expansion of JW n (z, A): 



W n (z, A) = y/W n M + £ (47) 

P >i A 

as well as the asymptotic expansions of the limit loci of zeros £^(A), I = 0, ±1, ±2, 
of the FS ip k (z, A): 



C (fe) (A) 



r (fc) 

Z^n>0 Z^ o =0 \(m + l)p+o ' t — 



■jp>0 Z~iq=0 \(m + \)p+q 

A.k) 

l^ P >0 l^q=0 A (2m+3)p+ 9 "T 



1 ^2m+2 ^i;2 P +l,q _ 9 , 1 

2m + 3 Z^p>0 Z^g=0 \(2m+3)v + q 1 11 ~ ' L 

^ A 2* 



Having done properly the semiclassical expansions of the respective quantities in 
order to be as close as possible to the previous formulations of Theorem 2a and 2b, 
we can come back to the previous form of A-dependence by a substitution in all the 

n _|_2 

above formulae A "a - back by A itself. Then we can formulate the following theorems 
analogous to Theorem 2a and 2b. 
Theorem 3a 

In the non-critical case and in the regular limit A — > oo zeros C/gr(A), |A| = [|A|] + 
A, / = 0, ±1, ±2, q = 1, 2, r = 0, ±1, ±2, of ij) h (z, |A|e i/3 j, < \(3\ < I are 
distributed on C cut uniquely along the corresponding exceptional SL's according to the 
formulae: 



L« W) - ^(y. A)) iy = ± ( 9 [|A|] + r - \) £ (49) 



'1 /77T^ — 1 „ , ,s\ , . / 1\ «7T 

where Ki((i k J r (\)) is a contour which starts and ends at Cz^r(A) rounding the turning 
point zi anticlockwise. 

The two lowest terms of the semiclassical expansions of zeros Ci^rW 9™ en by (48) 
are the following: 

L , k) , \ ^W n , a (y)dy = y/W n ,Mdy = ±qi*e-* 

<SU>W = Mr-*-i> , (50) 

V VV n,a{(,l, qr ;0,o) 

and are calculated from (51) according to the formula (68) o/App.2. 

The signs ± above are to be chosen to agree with signs of 3 (e tl3 /* \JW n (y, tydyj 
on the corresponding exceptional lines. 

Suppose now that there is an inner SL linking the roots Zk , z_k of W n (z, A) while 
the others are absent (see Fig.5b). In the limit case A — > oo it is possible only for 
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argA 7^ but with argA ~ |A|~«+2 (see App.3). Then for the assumed arrangement of 
SL's we have the following: 
Theorem 3b 

In the critical case, when there is inner SL linking z^ with z_k , and in the regular 
limit A -»• oo zeros C£(», |A| = [|A|] + A, I = 0,±1,±2,..., q = 1,2,..., r = 
0, ±1, ±2, ofipk(z, \X\e tl3 ), < \(3\ < -, k ^ k , —k , are distributed on C cut uniquely 
along the corresponding exceptional SL 's according to the formulae: 



L nm) (iv^> - a >) * = ± (*u + r - 



1 \ in 

a) X 



(51) 



where Ki(([ k J r (\)) is a contour which starts and ends at Ci'qrity rounding the turning 
point z\ anticlockwise. 
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Fig. 5a The regular limit A — > oo of zeros of ipk{z) for the potential 
W n {z) + 1 (n - odd). The non-critical case. The bold lines are exceptional SL's 

The two lowest terms of the semiclassical expansions of zeros Qor(^) 9^ ven by (48) 
are the following: 



A*) 

W,<jr;0,0 



W n , a (y)dy = ±qme 



-0 



Co,i(A) = ±(r - qA - -) 



me 



(52) 
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and are calculated from (51) according to the formula (68) o/App.2. 

The signs ± above are to be chosen to agree with signs of 3 (e zl3 f^ \J\V n (y, \)dy 
on the corresponding exceptional lines. 

In the cases I = k , —k the numbers q are bounded, i.e. q < |ifc |/7r with the minus 
sign in (52) chosen by assumption. 

Additionally in a regular limit Afc 0iS — > oo ; i.e. with R fixed, where \k 0:S — —^j-^i^ = 

[^k ,s] + Afc 0)S (-R), s = 0,1,2,..., there are two infinite sequences of zeros C-kl%-> 1 = 
1, 2, r = 0, ±1, ±2, of 'ipk ( z , -M distributed along the two SL's of the sector S-k 
according to the following rules: 

- - +r -\ + i)it, + ik, k Zk " dy ± 

JL- In 2 cos (ftr + ±9 Z h4 y) ± (d&A (53) 

■urat/i i/ie following lowest coefficients of the corresponding semiclassical expansion of 
Ak )± . 

*»— ko,qr' 



c (fe )± 



j^-k 
Jz-kn 



°'° \jW n , a (y)dy = -qin 



Kg 

C3&*,(*) = - (r - 9 A,,„(fl) - i + I T i ln2 cos(^)) * (54) 

V W n,a{^-k ,qr;0,o) 



where the plus sign corresponds to the SL being the upper boundary of S-k while the 
minus one to its lower boundary. 

Obviously the second part of this theorem applies also to the FS ip-k (z,\) with 
appropriate changes. 



7 High energy behaviour of zeros for quantized en- 
ergy 

By a quantization of energy we mean an identification of some two chosen funda- 
mental solutions. However such a choice is in fact quite limited particularly when the 
high energy limit, |A| — > oo, is taken. 

First, in general, one can not identify two neighboured FS's since they are always 
linearly independent. 

Secondly identifying in the high energy limit case any two not neighboured FS's 
we demand, as we will see below, the sectors the fundamental solutions are defined 
in to be connected by an inner SL. This means that we can consider only the stan- 
dard SG cases, i.e. when these inner SL's are parallel to the real axis of the C cut - 
plane. Therefore we can put the energy quantization conditions only on the pairs 
i>a,k{ z -> A), ip a -k(z, A), k = ±1, ±2, for the potential W n ^ a {z) and for analogous 
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pairs for the potential W n (z, A) + 1. But as it follows from each of the formulae (83) 
and (84) of Appendix 2 the corresponding quantization conditions for the potential 
W n ,a(z) when matching the pair ip a ,k(z, A), ip a -k(z, A) read: 



1 + exp 



a k \JW nja (y) + Z k (y, A) dy 



or 



°k^s\/W n!a (y) + Z k (y, A s ) ) dy = -(s + -)iri, s = 0, 1, 2, 



(55) 



(56) 



what proves that in the limit A s — > oo the SL emerging from z k ends at i.e. this 
is the inner SL. 
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Fig. 5b The regular limit A — ► oo of zeros of ipk{z) for the potential 
W n (z) + 1 (n - odd). The critical case. The bold lines are exceptional SL's 

The eq.(56) thus proves that for the potential W ntCl (z) when the energy is quantized 
there are inner SL's for each k = ±1, ±2, but only for one value of k, i.e. this 
which satisfies (56), the energy is quantized and the corresponding FS's coincide, i.e. 
ipa,k( z : A s ) = i(— l-) s ipa-k( z , K) as h follows from (83) or (84). However in the case of 
the potential W n (z, A) + 1 the SL satisfying eq.(56) can be the unique inner one for the 
corresponding SG while the remaining SL's of this graph emerge to infinity. 

Let now eq.(56) be satisfied for k = k . Then compairing (56) with eq.(34) we see 
that R = | for this ko and therefore we can not apply the result given by (34) to this 
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case. Nevertheless we can use the fact that now ^.^(z, A) coincides with ip a ^ (z, A) 
up to a constant and since it can not vanish in its sector S-k so does ip at k (z, A), i.e. 
there are no longer roots of ip aj k {z, A) along the boundaries of the sector S-k - 

Nevertheless there are still such roots of ip a ,k{z, A) in the corresponding sectors S-k 
for k 7^ k . 

Let us call the fundamental solution ip a> k (z, A s ) quantized if the quantization condi- 
tions (56) is satisfied just for the number ko. Of course it is satisfied also for the number 
— fco so the solution ip a _k (z, A s ) is also quantized and both the solutions coincide up 
to a constant. 

Therefore if we identify the quantized solutions 4> a ,k (z, \ s ) and ip a _k (z, X s ) then 
an exceptional set for them is just a section Ur=i ^ o nUr=i ^-k while for the remaining 
FS's their exceptional sets are kept unchanged. 

We thus come up to the following conclusion: 

A Imz 




Rez 



Fig. 5c The singular (R = |) limit A r — > oo of zeros of ipk (z) for the potential 
W nta (z) (n - odd). The critical case. The bold lines are exceptional SL's 

Corollary la 

In the case of the potential W n ^ a (z) when the singular high energy limit \ s —>■ oo 
(or equivalently s —>■ oo) are considered roots of both the quantized FS's and the not 
quantized ones are distributed on the C cut -plane uniquely on the exceptional SL's cor- 
responding to these solutions and the corresponding formulae of Theorem 2b, where 
A and A should be substituted by X s and A s respectievly, are valid for these distribu- 
tions excluding the formulae (34) and (35) which are nor longer valid for the quantized 
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solutions (see Fig. 5c). 

Exactly the same notes as above can be done with respect to the quantized solutions 
ipk {.z,\ s ) = i(— l) s ip-k (z, A s ) corresponding to the potential W n (z, A) + 1. They are 
now both deprived of zeros lying on the infinite SL's emerging from the turning points 
Zfc (A) and z_k (\) while keeping their zeros distributed on the inner SL linking these 
two turning points and on their remaining exceptional SL's. 

It is therefore clear how the corresponding conclusion for the quantized high energy 
limit for the potential W n (z, A) + 1 should sound. 

Corollary lb 

In the singular high energy limit X s — > oo roots of FS's for the potential W n (z, A) + 1 
are distributed uniquely on exceptional lines corresponding to these solutions according 
to the formulae of Theorems 3b, where A and A should be substituted by A s and A s 
respectievly, except the quantized solutions ipk {z, A s ) = i(— l) s, 0_fc o (z, A s ) for which the 
formulae (53) and (54) are nor longer valid (see Fig.5d). 




Fig.5d The singular (R = limit A r — > oo of zeros of ipk ( z ) f° r the potential 
W n (z) + 1 (n - odd). The critical case. The bold lines are exceptional SL's 

8 Summary and discussion 

In this paper we have shown that the high energy limit distributions of zeros of 
appropriately scaled fundamental solutions for polynomial potentials can be described 
completely both for the quantized and non-quantized cases of these solutions. The 
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quantized cases were considered earlier by Eremenko et al [9] where the authors pointed 
out that the ESL's are the loci of zeros of FS's. We have completed their observation in 
this case giving a detailed description of positions of these zeros on the ESL's. However 
we have considered the zeros distribution problem of FS's in general showing that loci 
of zeros of FS's on their ESL's is their common property which has to be completed 
by Theorems 2b-3b which find additional zeros of FS's outside the ESL's for the 
crititical cases of SG's. 

However, to get stable patterns for zeros loci distributions of FS's we have been 
forced to consider regular asymptotic high energy limit of these loci getting as a result 
island pictures for these distributions with the islands numbered by q in the correspond- 
ing theorems. Different regular limits (controlled by the A and R parameters) have 
lead do different distributions of zeros inside each island with the latter distributions 
being controlled by the r-parameter in Theorems 2-3. 

It is worth to note also that taking these regular limits stabilizes the limit zeros 
distributions of FS's in the same way as the rescaling z- variable by energy E in the 
initial polynomial potential P n (z) stabilizes its zeros distribution in the limit E — > oo 
reducing it to the zeros loci of the polynomial {—iaz) n — 1. 
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Appendix 1 

We shall show here that the integrals (6) are all pure imaginary. We simply 
calculate them. To this goal it is necessary to consider W n>a (z) on the cut z-plane. We 
make all these cuts from each turning point z^ parallelly to the real axis and running 
to ^tz = — oo if zj, lies to the left from the imaginary axis or is placed on it and running 
to 3lz = +00 in the opposite case. Then argument of z — z k will be taken from the 
interval (— it, tt) for the first group of turning points taking the values — tt, it below and 
above the corresponding cuts respectively (see Fig. la). For the secod group of z k their 
arguments are taken from the interval (0, 2ir) with the values 0, 2tt above and below 
the cuts respectively. 

Now, making in the integral (6) the following change of variable z — > £ = (— iaz) n 
we are to calculate the integrals 

^W n , a (z)dz = f ^W n , a (z)dz = -— [ (f - ifH^di (57) 

where the integration contour C k is an arc of a unit radius starting from the point 
z k and ending at z — z^k, running clockwise and avoiding all the met turning ponts 
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Zk-i, ■ z-k+i (see Fig. la) while C& is an image of Ck on the £-Riemann surface defined 
by the above change of variable (see Fig. 6). 
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Fig. 6 2n cut £ -planes of Riemann surface of (£ — lja^ - ^. The letters ai, d n 
denote boundaries of cuts. Boundaries with the same letters are glued 

It is easy to note that Ck is a collection of 2k or 2k — 1 unit radius circles on this 
£-Riemann surface depending on a. The first value corresponds to a = 1 and the 
second to a = e". The first circle starts at the point £ = 1 on the sheet on which 
arg£ = 2kir. The latter argument is just the one of £fc = {—iazk) n . In Fig. 6 this is the 
fc-th sheet from above with a^+i as the lower boundary of the cut between and 1 on 
this sheet. Each next circle is an image of successive arcs of the contour Ck- The last 
circle corresponds to the unit radius arc between the points Z-k+i and Z-k and it ends 
at unity on the n — k + 1-th sheet. 

It is easy to note that Ck is a collection of 2k or 2k — 1 unit radius circles on this 
£-Riemann surface depending on a. The first value corresponds to a = 1 and the 
second to a = e~~ . The first circle starts at the point £ = 1 on the sheet on which 
arg£ = 2kir. The latter argument is just the one of = (—iazk) n . In Fig. 6 this is the 
fc-th sheet from above with a^+i as the lower boundary of the cut between and 1 on 
this sheet. Each next circle is an image of successive arcs of the contour Ck- The last 
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circle corresponds to the unit radius arc between the points z_k+i and z_k and it ends 
at unity on the n — k + 1-th sheet. 

As it follows from the Fig.6 the integration along Ck on the £-Riemann surface can 
be deformed to go along the cuts between the points £ = and £ = 1 on the respective 
sheets of the surface. Then the integrations along the a^i and a„_ fe+1 cuts survive only 
(the remaining ones mutually cancel) so that their contributions to the last integral in 
(57) (call it ifc) are followings 



h = -— i (e-i)V^ 

ian Jc h 



-^/o(! -0 2 C » ^ ( -e l 2 - +e l 2„ ) for a = \ 



n n 



2i 



B(l±) for a = l, k = l,...,k<* 
sin (|, I) /or a = e-v, fc = l,...,fc< 



(58) 



where B(x,y) is the betha function. 



Appendix 2 



Here we prove the Lemma of sec. 5, Theorems 2a-2b of sec. 5 and Lemma' of 
sec. 6. 

Proof of Lemma 

To prove this lemma let us note that uj a (z) is holomorphic in C cut (e) vanishing there 
as z~2 n ~ 2 for z — > oo. Therefore the integrals Ik(z) = f lk ( z ) \w a (£)d£\i k — 1, ...n + 2, 
are well defined in the closure C cut (e) and are bounded there. Hence C a;e does exist and 
is finite in C cut (e) and therefore can be taken the same for each D^^. The estimation 



(15) follows then directly from (13) since 



1 _ e -2a k \w n , a (yi,y i+ i) < 2 in this formula for 
each % if all the integrations are performed on canonical paths jk( z ) what is allways 
possible by definition of Dk >e . 
Proof of Lemma' 

As we have mentioned in sec. 6 the turning points Zk(X) of W n (z, A) tend to corre- 
sponding roots z k of W nya (z) in the limit A — > oo. Therefore taking A > sufficiently 
large we can find all the roots Zk(X) in the corresponding circles \z — Zk\ < e for |A| > Ao- 
Therefore as in the case oiu> a (z) also w(z, A) corresponding to the potential W n (z, A) is 
holomorphic in the closure C cu t(e) and bounded there since it vanishes there as z~z n ~ 2 
for z — > oo. Hence, as previously, the integrals Ik(z,X) = J 7fe ( 2 ) A)d£|, k = 

l,...n + 2, are well defined in the closure C cut (e) and are bounded there by C e (A) 
defined by 

C e (A) = liminf 4(^,A) (59) 

7 fc (z),zeC cui (e), fe=l,...,n+2 
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However C e (A) is a continuous function of KA and SA for |A| > Ao- Therefore since 
also liniA-^oo C e (A) = C a - e in a closure of C cut (e), C e (A) is bounded by some constant 
C e (Ao) with the property that C e (Ao) — > C a;e for Ao — > oo. The latter means that for 
Ao sufficiently large we can choose a constant C e (Ao) to be independent of Ao and to 
differ from C a . t by e, i.e. C e (A ) < C e = C a . t + e. 

Now the estimation (38) follows directly from (12) written for the potential W n (z, A) + 

1. 

Proof of Theorem 2a 

According to Lemma, for every k, ip a ,k(z) can be represented in D k , e for |A| > 
A > C k)6 as 

lM*) = ^(*)^ 7 >1 (60) 

In fact the above estimation of ip a ,k{z) in -Dfc,e can be extended to any finite order 
in A -1 . It will be convenient to include into this estimation all orders of A -1 i.e. to 
represent ip a ,k{z) by its full asymptotic form (16). Therefore we will use further the 
following asymptotic representations for ip a ^{. z ): 

Mz) ~ r a %(z) = W~kz)e kX ^ ^^ dy x t°( z , A) = 

W~! {z)e kX ^ VWn, a (y)dy+C k z k ( y ,x)d y ^ 

Dispite the fact that the asymptotic semiclassical series Z k (y, A) is generally diver- 
gent we can use it in its full form having however in mind that it can be abbreviated 
in any moment at some finite power of A -1 to provide us with an estimation like (60) 
but valid then up to the order kept. 

To perform detailed calculations we shall assume the C cut -plane to be cut in the 
following way. 

For the root zq and z™ (if they are) the corresponding cuts coincide with the positive 
and negative parts of the imaginary axis respectively. So there are in fact two SL's, 
mutually parallel to themselves on C cut and lying on different sites of the cut. 

For the remaining roots the corresponding cuts are parallel to the real axis running 
to the left for the roots lying on the left from the imaginary axis and runnig to the 
right for the roots placed on the right from the imaginary axis so that the pattern of 
cuts is completely symmetric with respect to the imaginary axis. 

Arguments of differences z — z k , k — 0, ±1, ...,n/2, are always taken with respect 
to axis emerging from the roots z k parallelly to the real axis and keeping its direction 
(see Fig. la). For the right hand roots these arguments are taken from above the cuts. 
The arguments are positive if they are taken clockwise and negative in the opposite 
cases. 

With these conventions a sign of & (A J z z k JWZMdy) , | arg A| < ^ is therefore 
always positive for the sectors So, S ± n+3 and S ± n+2. This sign is always negative for 
the sectors Sn+i and S±«. For the remaining sectors S k , k = ±1,±2, this sign is 
negative above the cuts crossing these sectors and positive below them, i.e. it changes 
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each time these cuts are crossed (see Fig. la). Note that this switch of the signs does 
not occure on the other cuts. In the same way behave all terms of the expansion Zk in 
(61). 

_ i 

Also the factor W n ,a{z) in (61) crossing each cut changes its phase by =p. The 
upper sign has to be taken when the cut on the positive imaginary axis is crossed from 
its right or the left hand cuts are crossed from above while the lower has to be taken 
for the negative imaginary axis cut and the r.h. cuts crossed in the same directions 
(i.e. to the left and down respectively). The signs are changed to opposite if the cuts 
are crossed in the opposite directions (see Fig. la-lb). 

Therefore we assume the form (61) of the solution ip a ,k{z) to be valid in the negative 
parts of the sectors S&, k = ±1,±2, .... and in the sectors Sn+i and S,n+2, so that 

2 - 1 - 2 

a k = +1 there. 

The form (61) is also assumed to be valid for the sectors So, S ± n+i and S±r but 

then <7fe = — 1 for the corresponding solutions. 

The asymptotic forms (61) of the FS's when the latter are continued analytically 

across the C cu t-plane along canonical paths while keeping their validity change therefore 

_i 

appropriately by acquiring additional phases by the factor W n ,a (z) or by switching their 



signs by the exponentials A J z 2 fc \JW n ^ a (y)dy + f£ Ck Zk(y, X)dy on the cuts met. 

We shall consider below in details the analytic continuation of the exact solutions 
?pa,k(z) and their asymptotic expansions (61) as well along canonical paths. We will 
get first exact results of such continuations to come over next to their asymptotic 
forms and to get finally their limit for A — > oo. We shall consider them one by one 
beginning with k — 0, 1,2, i.e. going down the corresponding SG's, 

but neglecting detailed calculations for the solutions for which the corresponding results 
can be obtained be symmetry arguments. 

k — case 

This case corresponds to an even-n and a = e~~ so that the corresponding SG 
with exceptional SL's for this case is shown in Fig.ld-le. 

We shall continue ip a ,o(z) to vicinities Vk, e of its exceptional SL's (ESL) expressing 
it simply as a linear combinations of FS's defined in sectors closest to the particular 
ESL considered at the moment. By symmetry arguments we can do it only for the left 
ESL's of Fig.ld-le. 

Using first the exact forms (11) of the solutions we can get coefficients of the 
coresponding linear combinations also in their limit forms up to any order of A -1 
by substituting the exact expressions by their asymptotics (61). 

Therefore continuing ip a ,o(z) to the vicinity of the ESL emerging from the root 
z k , k = 1, f we get: 



Xk^k+l 
Xk^k+l 
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Wn,a{z)e Jz o V / -A f Jw^fidy , 



-«e Jz fc v Xo-kXk+i(*» (62) 



where x^j = hm^^ Xi(^) = hm^^. Xj( z ) = Xj^i (see, for example [4], ref.2). 

From (62) we get the following exact condition for loci (j^ of zeros of ip a fi{z) on 
the considered ESL: 

I ( 1 Xo^fc+iX^C^i.A) 

where m is positive integer and the choice of the "+"-sign is determined by the posi- 
tiveness of the l.h.s. of (63) on the ESL considered. 

Now we would like to come over in (63) to the limit A — > oo. However if we want 
to use the results of Lemma we need C^ m to be kept inside the domain D ky€ , i.e. their 
distances to z k in the limit A — > oo should be greater than e. It then follows from (63) 
that for large A also m should be large, i.e. 



mir 

— — > hmsup 
A 



<7T 



z k +ee 



%4> 



W n , a (y)dy 



4(e) (64) 



This condition means that we have to take into account only those zeros (j^ which 
are sufficiently far from z k . Therefore for those which can fall onto z k we have to take 
m sufficiently large, i.e. m > m Q = |A|ifc(e)/7r. 

For the remaining ones, i.e. those which are to approach finite non-zero distances 
to Zk it is necessary to make m growing linearly with A, i.e. m = q[\X\] + r with 
q = 1,2,..., r = 0,±1,±2,..., and [|A|] being the step function of |A|, i.e. |A| = 
[|A|] + A, 0< A< 1. 

Let us note further that the latter case will cover also the previous one when q = 
and r > mo- In all the formulae below we shall assume the cases of zeros falling down 
onto Zk to be taken into account just in this way. 

Having this in mind we can come in (63) to the limit [|A|] — > oo to get for ~ 

CfcJr(A): 

f A 0) (A) , / 1\ .V 1 f 

I" >WWA' = + («+'--l)T + 2A / k>(C(a)) ^fe.A)* (65) 

where KkiCklri^)) * s a con tour shown in Fig.4a which starts and ends at the point 
C|°] r (A). This contour is not closed since it starts and ends at different sites of the cut 
emerging from Zk- 

In fact (65) is an implicit condition for the asymptotic expansion (j^qrW °f zeros 
Cfc°„ as A — > oo. Therefore we can look for the following form of semiclassical expansion 
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for OA): 



OA) = E ^C(A) (66) 



p>0 

with C^ r0 (A) = z k . 

Let us now note that the condition (65) can be written uniformly as: 

/«<>„ G^vw - ^-to. *)) * = + (riiAii + - j) x < 67 > 

We can now calculate the asymptotic series (66) in the limit [|A|] — > oo with fixed 
A from the following formula: 



+ 



2 E^ (^7wv(y)-^a,o(y,A) 



E/ ^ p Cfc,qr-p(A) 



(^[|A|]+r-i)^, g>0 (68) 



and with a similar formula for q = obtained from the last one where Cfc°g r o(A) is 

substituted by Cfc°o r i(A)/A and p > 1. 

The limit []A|] — > oo with fixed A (and argA = (3, \(3\ < n/n) considered above is 
of course regular. 

Therefore for q > the zero term C^ qr0 {A.) of such an expansion is given by the 
equation: 

L , A „ yW n , a (y)dy = [ <iUA) y/W^Mdy = qiire^ (69) 

which shows that this term is independent of A. 
The next term is given explicitely by: 



and so on. 

For q = we get correspondingly 



and 



/ V / MVWy = (r- i ) y , r > m = |A|4(e)/7r (71) 
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k = ^ case 

This case corresponds to n odd and to the pattern of ESL's shown in Fig. la. We 
have to continue ib a ^ n+z (z) to vicinities of ESL's emerging from the root Zq and from 

the roots z k , k = ±1, ±2, i 23 ^. In the first case we have to express ib a ™+s (z) 
linearly by the solutions ib a ™+z (z) and tp a _i(z). For the remaining ESL's ib n+z (z) 
is continued in the same way as in the previous case. We have therefore: 

2 ^-2+2(2) +ie J *o v — ? a-i(^) (73) 

for the first ESL and 

Xk^k+1 

ie -xQ V^^fir ^^ dy ^^^ k+1 (z) (74) 

Xk^k+1 

for the ESL's emerging from the roots z k) k = 1, 2, while for the ES's emerging 
from the roots with the opposite sign of k we get: 

X-k-+-k-l 

. e -A/;- v^-^^a/-- v^)^ (75) 

X-fc->-fc-i 

Arguing in exactly the same way as in the previous case from (73) and (75) for 
zeros of ip a n+3 (z) on the ESL's emerging from the roots z± k , k = 0, 1, 2, 21 =^ we get 
the corresponding conditions: 

^<df>>, (r/^w ~ as*** x) ) dy=± ( m +r -i)f (T6) 

= ^ case 

This case corresponds to the even n and a = 1, see Fig. lb. We have to continue 
ip a n±i(z) to vicinities of its the ESL's emerging from the roots z± k , k = 0, 1, 
and from za. 

2 

However, the corresponding formulae for the positions of zeros on the ESL's emerg- 
ing from all the roots mentioned above except the last one coincide for the obvious 
reasons exactly with the respective formulae (76) while for the ESL's emerging from 
the last root we have 

r.z n J 

ie -*Sj 0*,.G0* ( Xaj p_^ a> .(*) - x^^(z)) (77) 
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where z is assumed to lie on the l.h.s. of the cut emerging from z«. 

(2+2) 

Therefore for zeros ( R 2 qT of ip a n+2 (z) lying in vicinity of the ESL emerging from 
zn we get: 

L«P>m («[W1 + r-i) * (78) 

where ± corresponds to the positions of (n. qr on the left or right from the cut re- 

spectively. Of course we always get the same values for £n 2 qr irrespectively of the site 
chosen. 

k = 1, 2, | cases 

We shall consider all the enumerated cases jointly since the corresponding roots are 
all on the left from the imaginary axis independently of the parity of n. 

We have to distinguish two cases corresponding to arg A = (3 is equal to or differs 
from it. In the second case (the corresponding SG is shown in Fig. 2) the procedure 
of analytical continuations of ip a ,k( z ) along canonical paths to vicinities of all its the 
ESL's is completely analogous to the previous cases since all the sectors having these 
ESL's as their boundaries are available from the sector S k along canonical paths. The 
rules of making these continuations are now clear and their results are gathered by the 
formulae of the type (67), (76) or (78) so we can focuse our attention on the real case 
of A. 

In the case (3 = each inner SL emerging from the root z k) k = 1,2,...,| is 
exceptional for the corresponding solution ipa t k(z), see Fig.la-le, together with the two 
ESL's which close the sector S- k which is not available directly on any canonical paths 
from the sector S k . 

On the other hand positions of other ESL's corresponding to ip a ,k(z) m the consid- 
ered case is analogous to their positions for (3 ^ and their vicinities can be approached 
by continuations of ip a ,k( z ) along canonical paths, i.e. positions of zeros of the con- 
sidered solution on these ELS's can be established exactly in the same way as for the 
cases investigated so far. Therefore for positions of these zeros on their ESL's we can 
invoke again the respective formulae (67), (76) and (78). 

Consider now positions of possible zeros of ip a ,k( z ) on the inner line emerging from 
the root z k (it ends at z^ k ) and on the remaining two lines emerging from Z- k and 
runnig to infinity. For the inner line we can use any pair of linear independent FS's 
for both of which this line is not exceptional. These can be for example the solutions 
ip a -k+i(z) and ip a ,-k-i( z )- They communicate canonically with the solution ip a ,k( z ) s° 
we get: 

X-fc+l->- k-1 
A r _fc_1 JW n , a (y)dy Xk^-k+1 i i x 

X-k+l^-k-1 
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(79) 



From (79) we get for the positions of zeros of if) ai k(z) on its inner ESL (zk, Z-k)' 

- k 2 ^ A) ) * = " ( m " i) t (80) 



It follows from (69) that q is limited now by the integral Ik = f* k k yW n , a (y)dy, i.e. 
q < -7fc/(i7r). 

In order to look for possible zeros of ip a ,k{z) on its ESL's emerging from we have 
to continue it further to vicinities of these lines using the solution %p a _k{z) which also 
comunicates canonically with both the solutions ip a -k+i{z) and ^ a -k-i{z). Choosing 
one of them we can express it as linear combinations of the others to get: 



ip a - k +i(z) = 



6 k X-k+l^-k-lVa,-k{Z) 



il>«,-k-i(z) (81) 



and 



ip a -k-i(z) = 

k+1% 

-a r- k+1 Jw n , a (y)dy-^ f z - k - 1 Jw n , a {y)dy 



£- k -' ^w n , a {y)dy , , 

e z - k X-k-l^-k+liPa-k{z) 



ip a ,- k+ i(z) (82) 
respectively. 

Now substituting subsequently both the above formulae to (79) we get the formulae 
realizing continuations of ^ a , k {z) close to the respective ELS's, i.e. we get in this way: 

, / \ ~ A f"~ k ^w n , a {y) d y Xk-*-k-i , i ( \ 

Wa,k{z) = -e Jz k v (x-fc+i— fc-iVa -k{z)- 

X-k+i->-k-i 

\ Xk-^-k-l) J 

and 

, / \ , A f"~ k \f w ^Ay)dy Xk^-k+i , I / \ 

Va,k{z) = +e Jzk V (X-fe+l— fc-l^a ,-k{Z)- 

X-k+l^-k-1 

V Xk^-k+iJ J 

The first of these formulae is suitable for analysing zeros of tp a! k(z) in a vicinity of 
the SL being the lower boundary of the sector S-k while the second in vicinity of the 
SL being its upper boundary. 
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The first (83) of these formulae provide us with the following condition for zeros 
C-l~m °f i'a^iz) in the limit A — > oo: 



/ C ; m y/W nta {y)dy = - (m - i) ^ + -L (A ^Ji)dy + ^ Z a>Jfc dy 

In 2 cos ^ (A / sJW n , a {y)dy + / Z a>k dy) -—( Z a , k dy (85) 

where ^ and K_ k {C < ^l~ m ) are contours shown in Fig. 

Putting now A rfe JW ny0l {y)dy = 2{r k + i?)?7r, \R\ < \ with a non-negative 



integer, A rfc = [A r J + A rk (R), < A rfc < 1 and m = q[\ rk ] + r we get from (85): 

v^fei* = - («w + r - 1 + f ) ^ + ix; 1 ^ - 

-L In 2 cos (flx + I» Z, A ) - -L (86 ) 

where the limit taken is regular with fixed R. 

Therefore for the first two terms of the semi classical expansion of C-l~qr we g e ^ : 

/ k '" ri) Jw n , a (y)dy = -qm 
C%i(R) = - (r - qK k {R) - - + f + \ ln2cos(i?7r)) t= (87) 

In exactly the same way we get the condition for zeros (-I'm distributed along the 
upper SL emerging from the turning point Z- k in the limit A — > oo: 

( .(k)+ j . 

IT l/^O* = - (•» - j) T + 4A 1 + ^ + 

i- In 2 cos i» (A / <jw v (y)dy + / Z^dy) + -L / Z„, fc A/ (88) 

and with the same meaning of notations used we get: 

(fc) + 

J- In 2 cos (flx + i» ^ Z a ,iy) + J- Z<a d» (89) 

with the following first coefficients of the corresponding semiclassical expansion of 

Ak)+ . 
t>— k,qr' 

fc,qrO 



\fW n , a (y)dy = -qm 

eiUR) = ~(r- qK(R) + ln2cos(^)) ™ (90) 
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k = ! ^ case 

This is the last distinguished case which has to be treated. It corresponds to an 
odd n of Fig. la and to the solution il> a ^+i (z). We should continue the latter solution 

to its ESL's emerging from the roots Zk, k = 1, t± y- and from zq. Continuing to any 
ESL of the first group we get: 

A f fe JW n , a {y)dy-\ P*" 1 Jw n , a {y)dy i A P fc ^"."(j/)^ 

e X2±i-itW-i(2) = Wn,a{z)e -r- x 



|e J ^ v Xn±i^ fc _iXfe(^) X2±i^fcXfc-i(^) 1(91) 

and hence for the respective distributions of zeros ( k q 2 r (A): 

(s^- tx z «* (i >' A) ) * = ( m] +r -M (92) 

For the ESL emerging from zq we have: 



le 



(xn±l^n±3Va,2±2(^) -X"±l-»-»±3V'a-H±2(^)) = 

i A P° \/W„,a(y)(iy / r * /ttt — rr , 
lWn,a{z)e T- le Jz o v y^ 1 : "+3Yn+3(z) + 



«e Jz o v Xn+i-^-ii+aX-ii+i l- 2 ) ) (93) 



2 2 2 

and hence: 

L&> m G^ra - A >) * = («n A| i + r - 1) x < 94 > 

for the respective distributions of zeros on the considered ESL. 



Appendix 3 

We show here that for the potential W n (z,X) in the limit A — > oo if there is a 
unique inner SL of the SG for this potential between the turning points Zk (X) and 

2 

Z-k (X) then argA 7^ and in general case argA ~ |A| n + 2 . To see this consider in the 
limit A — > 00 the integral A /^"(^ \JWrJjJ, X)dy. With the accuracy O (|A|~"+2^ we 
have: 



A / yJW n (y, \)dy = A / yjw a>n (y)dy 



— AO-A—ia) / ; — — (95) 

2n J ^ JwZmx^ 
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where we have used a relation Zk(X) = Zk — + 0(X~^) valid for any 

turning point in this limit, with Zk being a turning point of W atTl (z). 
According to our assumption we have further (A = |A|e l/3 ): 

ft U [ Z ~ koW y/w n (y, Xjdy) = - J-» K_iH«) B_2 f "° ^==7 777^(96) 
V \W J 2n \ J **o yJW a>n (y)dy) |A|^ 

and hence finally: 

P = -7T^^ ( K l(-^) n " 2 I*" 1 = ) A- (97) 

where 4 = //-*o yjw n , a (y)dy. 
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